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ABSTRACT 

The  general  equations  of  motion  for  a  rigid  body  are 
derived  in  cylindrical  coordinates  by  Lagrangian  dynamics 
and  used  to  model  the  motion  of  an  airplane  in  a  steady  spin. 
After  simplification,  the  equations  are  cast  into  a  form 
utilizing  conventional  aerodynamic  derivatives  along  with 
other  derivatives  which  may  be  significant  in  spins.   An 
iterative  numerical  solution  procedure  is  outlined  which 
should  simplify  the  problem  of  solving  the  nonlinear  dif- 
ferential equations,  and  relationships  between  the  Euler 
Angles  used  in  the  equations  and  the  more  familiar  ordered 
set  of  pitch,  roll,  and  yaw  are  derived  to  permit  computer 
input  and  output  of  orientation  to  be  more  easily  visualized. 
The  inverse  problem  is  also  considered ,  and  equations  are 
derived  to  translate  spin  test  data  into  the  six  Lagrange 
coordinates  used  so  that  computer  calculations  may  be  checked 
against  test  results  and  actual  aerodynamic  forces  can  be 
compared  with  computed  values.   The  paper  concludes  with  a 
discussion  of  the  implications  of  this  spin  analysis  along 
with  applications  and  extensions. 
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TABLE  OF  SYMBOLS 

The  definitions  of  the  symbols  used  in  this  paper  are  as 

follows : 

CQ  Center  of  gravity 

C      Mean  aerodynamic  chord ,  feet 

Cfc  C*   C/?  Aerodynamic  force  coefficient  along  the 

*  '*  subscripted  axis,  dimensionless 

^-f\\  (-M/j  Cm-  Aerodynamic  moment  coefficient  about  the 

subscripted  axis,  dimensionless 

F^j  Fyj^%  Aerodynamic  force  along  the  subscripted 

axis,  lbs. 

Fop.  Lagrange  generalized  non-conservative 

force,  lbs. 

Fit>t  fkt  />       Lagrange  generalized  force  in  the  subscripted 
/:  El     Fi..  coordinates,  lb.  or  ft  -  lbs. 

2 

Q  Acceleration  of  gravity,  ft/sec 

Q      ^vjQt  Airplane  eg  acceleration  components  along 

*'  the  subscripted  axis,  ft/sec2 

1^  Xy  X$  Airplane  principal  moments  of  inertia, 

slug  -  ft^ 

M  Airplane  mass,  slugs 

Mv.  Mv/  M-      Moments  about  the  subscripted  axis,  ft  -  lbs. 

Xy    fj  -C 

^       Sequential  subscript,  integer 

p t    G       \T  Angular  rates  of  rotation  about  the  airplane 

X,  Y,  and  Z  principal  axes,  respectively 

*\.Y  Lagrange  generalized  coordinate 

ft  Spin  radius,  feet 

v->  Airplane  wing  area,  ft 

7"  Lagrange  kinetic  energy  term,  ft  -  lbs. 

"t  Time,  seconds 


At        Increment  of  time,  seconds 

/        Lagrange  potential  energy  term,  ft  -  lbs. 

Ifo        Airplane  eg  velocity  ft/sec 

ML/n.  ^U/yV^W*        Airplane  eg  relative  wind  velocity  component 
'  along  the  subscripted  axis,  ft/sec 

X    y  2  Airplane  principal  axis  system 

Xtj  Yt  £,  Reference  Cartesian  axis  system 

E9  Airplane  altitude,  feet 

29  Rate  of  ascent ,  ft/sec 

oC  Airplane  angle  of  attack  relative  to  principal 

axis  system,  degrees 

°0j         Direction  cosines  in  terms  of  Euler  angles 
(see  Table  I) 

fi  Airplane  yaw  angle,  degrees 

Y  Lagrange  generalized  coordinate  indicating 

direction  of  spin  radius  vector,  radians 

o  Spin  rate,  radians/sec 

4l  **"  *h  Aileron,  rudder  and  elevator  deflections, 

degrees 

Q/  $  <p  Airplane  principal  axis  system  orientation 

Euler  angles  with  respect  to  the  X-,,  Y-,  , 
Z,  system,  radians 

Q/rfajrh  Airplane  orientation  in  pitch,  roll,  and 

yaw  with  respect  to  the  X-,  ,  Y-,  ,  Z,  system, 
radians 

&PjTA*fy  ,  Pitch,  roll  and  yaw  rates,  radians/sec 

F~     £v  Ez  Incremental  angular  rotations  about  the  X, , 

1'   1  sYstem  a^es/  radians 

f  Air  density,  slugs/ft 

CJy,^-)//^       Angular  rates  of  rotation  about  the  air- 
plane principal  axes,  equivalent  to  p,  q, 
and  r,  radians/sec 


I.   INTRODUCTION 

A.   FOREWORD  AND  ACKNOWLEDGEMENTS 

The  intent  of  this  work  was  to  develop  an  alternative  to 
the  conventional  method  of  computing  spin  characteristics 
which  thus  far  has  not  been  uniformly  successful  in  pre- 
dicting airplane  spin  modes.   In  this  effort,  the  author 
has  tried  to  approach  the  problem  with  the  capabilities  of 
modern  computers  in  mind  and  devise  a  method  of  solving 
directly  for  spin  modes  rather  than  mathematically  tracking 
airplane  motion  through  all  intervening  phases  of  motion 
preceding  the  spin.   If  any  success  is  realized  through  this 
approach  as  outlined  herein,  it  will  be  due  in  no  small 
measure  to  the  tireless  efforts  of  Lieutenant  R.  L.  Champoux, 
USN,  whose  contribution  of  months  of  hard  work  and  extensive 
knowledge  of  computer  programming  were  instrumental  in  the 
evolution  of  this  approach.   Lieutenant  Champoux" s  separate 
work  on  the  computer  aspects  of  this  analysis  may  be  found 
in  Ref.  9.   The  author  is  also  indebted  to  Professors  L.  V. 
Schmidt,  T.  H.  Gawain,  and  M.  H.  Redlin,  of  the  U.  S.  Naval 
Postgraduate  School  whose  counsel  helped  shape  many  aspects 
of  this  work. 

B.   BACKGROUND 

Airplane  spin  or  autogyration  is  a  mode  of  motion 
exhibited  by  virtually  all  fixed  wing  aircraft.   It  is  at 
best  an  undesirable  flight  maneuver  and,  with  growing 
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frequency  in  new  airplane  designs,  a  terminal  maneuver  from 
which  recovery  cannot  be  effected.   Although  the  flight  con- 
ditions which  lead  to  spins  can  be  avoided,  pilots  of  military 
light  attack  airplanes  and  especially  fighter  airplanes  do 
not  have  that  option  if  pressed  into  flying  to  the  "edge  of 
the  airplane's  flight  envelope"  by  the  mission  or  a  tactical 
situation. 

The  need  to  utilize  the  edge  of  the  flight  envelope  in  a 
combat  role  is  readily  apparent  when  one  considers  that  it  is 
only  in  this  fringe  area  that  the  superiority  of  one  airplane- 
pilot  combination  over  another  can  be  realized.   However,  less 
obvious  is  the  extent  of  additional  exposure  necessary  for 
training  in  this  area  if  a  pilot  is  to  become  proficient  in 
a  particular  airplane  under  such  conditions.   Thus,  flight 
near  the  edge  of  the  envelope  must  be  more  than  an  occasional 
event  for  the  military  pilot.   Unfortunately,  loss  rates  due 
to  irrecoverable  spin  modes  in  advanced  tactical  jet  airplanes 
frequently  makes  effective  training  in  this  area  economically 
prohibitive.   Thus  the  margin  of  superiority  built  into  the 
airplane  at  great  expense  is  in  many  cases  lost. 

It  is  in  recognition  of  this  fact,  and  the  fact  that  spin 
predictions  based  on  conventional  analysis  techniques  have 
not  shown  much  success,  that  this  effort  was  undertaken.   It 
is  not  expected  that  this  analysis  of  airplane  spin  will  lead 
directly  to  solutions  of  the  spin  problem  itself,  but  it  is 
hoped  that  this  work  will  cast  the  problem  in  a  new  mathe- 
matical perspective,  from  which  additional  insight  into  the 


spin  problem  might  be  gained  towards  a  goal  of  eventually 
designing  out  irrecoverable  spin  modes. 

C.   GENERAL 

Historically  airplane  spin  studies  have  employed  forms  of 
the  same  equations  of  motion  which  are  used  in  the  small 
perturbation  approach  to  airplane  dynamics.   This  method  of 
mathematically  modeling  the  spin  appeared  promising  because 
the  equations  could  be  readily  integrated  by  numerical 
methods  to  produce  time-histories  of  characteristic  parameters 
which  could  then  be  compared  directly  with  flight  test  time- 
histories.   It  was  also  apparent  that  the  body  axis  coordinate 
system  used  in  the  classical  approach  provided  simpler 
equations  for  random  tumbling  motion  than  those  obtained 
using  an  inertial  system  of  coordinates.   Thus  the  post- 
stall  and  incipient  spin  gyrations  of  the  airplane,  which  in 
many  cases  are  of  equal  concern,  can  be  more  readily  tracked 
by  a  computer  in  a  body  axis  system  than  an  inertial  system. 
The  degree  of  success  realized  by  thus  extending  the  classi- 
cal approach  however  depends  to  a  great  extent  on  the 
judgment  and  intuition  of  the  investigator  in  leaving  just 
enough  coupling  in  the  equations  to  permit  satisfactory 
representation  of  the  motion  without  introducing  excessive 
non-linearities.   Frequently  this  process  degenerated  to  a 
trial  and  error  parametric  approach  without  much  insight 
being  gained  and  without  general  validity  being  obtained  for 
all  airplane  configurations  or  spin  modes. 
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On  the  other  hand,  the  treatment  of  the  spin  problem  in 
an  inertial  coordinate  system  appears  to  add  undue  complexity 
to  the  computer  tracking  problem,  and  indeed  for  random 
tumbling  motion  the  non-linearities  become  excessive.   However, 
if  the  steady,  or  quasi-steady   spin  is  to  be  studied,  an 
inertial  cylindrical  system  offers  a  number  of  distinct 

advantages,  among  which  are  the  ability  to  differentiate 

2 
between  motion  due  to  spinning  and  motion  due  to  tumbling  , 

and  a  sufficient  reduction  in  the  complexity  of  the  coupling 

3 

to  permit  computer  solution  of  the  exact  equations  .   Since 

the  intent  of  this  study  was  to  develop  a  generalized  method 
of  solving  for  steady  spin  modes,  the  equations  of  motion 
were  derived  in  an  inertial  coordinate  system.   The  details 
of  this  coordinate  system  are  shown  in  Figure  1,  and  discussed 
in  some  detail  in  Section  III. 


Spinning  motion  in  which  state  variables  have  a  mean  com- 
ponent and  an  oscillatory  component. 

2 

Random  motion  inconsistent  with  steady  or  quasi-steady  spin 

motion. 

3 

No  small  quantity  approximations  are  made  in  either  the 

derivation  or  simplification  of  the  equations  of  motion. 
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II.   SCOPE 

This  work  was  preceded  by  a  survey  of  current  airplane 
spin  literature  [Ref.  1]  to  determine  the  present  state-of- 
the-art  in  airplane  spin  analysis.   From  that  study  it  was 
evident  that  possibly  another  avenue  of  approach  to  the 
problem  ought  to  be  explored.   Thus,  rather  than  attacking 
the  spin  problem  on  the  broad  formidable  front  of  mathe- 
matically following  an  airplane  through  all  transient 
gyrations  from  stall  to  spin,  it  was  decided  to  limit  the 
scope  of  the  investigation  and  probe  directly  into  the  final 
steady  state  motion  of  the  spin  where  a  force  and  moment 
equilibrium  situation  might  yield  to  simpler  mathematical 
representation.   In  this  approach  no  restrictions  were 
imposed  on  the  form  of  the  airplane  aerodynamic  data  to  be 
used  in  this  analysis  other  than  requiring  that  it  be  ref- 
erenced to  the  airplane's  center  of  gravity  and  principal 
axes.   By  leaving  the  details  of  the  makeup  of  such  data 
somewhat  general,  the  analysis  is  freed  from  any  incumber- 
ances  of  classical  aerodynamic  derivatives  and  thus  can 
proceed  more  directly  to  the  derivation  of  the  general 
equations  of  motion.   Although  conventional  force  and  moment 
derivatives  are  suggested  for  use  in  initial  computer  pro- 
gramming, subsequent  work  may  suggest  a  more  appropriate  form 
of  aerodynamic  force  and  moment  data,  or  possibly  more 
representative  aerodynamic  derivatives  may  be  devised. 
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III.   DERIVATION  OF  THE  EQUATIONS  OF  xMOTION 

The  coordinate  system  used  in  the  derivation  of  the 
equations  of  motion  is  a  combination  of  an  inertial 
cylindrical  system  and  the  airplane  principal  axis  system. 
The  vertical  (z)  axis  of  the  cylindrical  system  corresponds 
to  the  central  axis  of  the  motion  or  "spin  axis".   The 
advantage  of  the  cylindrical  system  is  that  it  locates  the 
airplane  center  of  gravity  in  convenient  terms  of  an  alti- 
tude coordiante  (z  ) ,  spin  radius  coordinate  (R)  and  angular 
position  coordinate  (ft) ,  when  the  airplane  is  in  a  steady- 
state  spin.   The  orientation  of  the  airplane  about  its 
center  of  gravity  with  respect  to  the  inertial  cylindrical 
system  is  specified  in  terms  of  Euler  angles.   A  cartesian 
coordinate  system  is  fixed  at  the  eg  position  on  the  (R) 
vector  having  its  (x, )  axis  in  the  (-R)  direction  and  its 
(z, )  axis  parallel  to  and  in  the  opposite  direction  of  (z  ) . 
This  cartesian  system  provides  a  reference  for  the  orientation 
Euler  angles.   Thus  for  zero  values  of  the  Euler  Angles  (9), 
(V)  ,    and  {(j))  ,  the  airplane  would  be  upright,  wings  level, 
and  have  the  (X)  axis  of  its  principal   axis  system  pointed 
inward,  along  the  (R)  vector,  directly  at  the  spin  axis. 
The  details  of  this  hybrid  coordinate  system  are  shown  in 
Fig.  1   along  with  the  positive  directions  of  the  position 
coordinates.   A  general  schematic  of  Euler  angle  relation- 
ships is  contained  in  Fig.  2  ,  and  equations  for  the 
direction  cosines  of  the  airplane  axes  expressed  in  terms  of 
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INERTIAL  REF. 
Fig.  1.   The  Coordinate  System  Used  to  Model  the  Spin 


Angular  relation- 
ships between  the 
airplane  x,  yj   z  and 
the  x{,  y     zt      axis 
systems  (shown 
inverted  to 
facilitate 
visualization) 


X1Y1  plane 
*A  negative  value 
of  <fi   has  been 
depicted  in  this 
schematic 


(INTERSECTION  OF 
X,-Y,  AND  X-Y 
PLANES) 


Fig.  2.   A  General  Schematic  of  Euler  Angle  Relationships 


14 


Euler  angles  are  listed  in  table  I.   The  rationale  behind 
this  choice  of  coordinate  systems  is  that  it  more  simply 
(in  a  mathematical  sense)  represents  the  motion  being 
modeled.   Thus  even  for  oscillatory  spins  the  motion  would 
be  readily  recognizable  and  translatable  into  mean  values 
of  the  six  coordinates  (three  position  and  three  orientation) 
In  addition,  since  the  principal  objective  is  to  solve  for 
steady  spin  modes,  the  complexity  of  stall  and  departure 
motion  in  this  coordinate  system  has  no  impact  on  the  results 

4 
Table  I:   Direction  Cosines  in  Terms  of  Euler  Angles 


Cosines  of  Angles   between  X,  Y,  Z   and   Xx,  Ylt  Zx 

x 

Y 

Z 

*i 

au    =    COS<£  cos  \f> 

—  sin  <j>  sin  <p  cos  8 

«2i    =    —  sin  <j>  cos  \f> 

—  cos^  sin  4>  cos  S 

<*3!    =    sin  8  sin  4> 

Y, 

oj2    =   cos  <f>  sin  \f> 

+  sin  <f>  cos  4>  cos  8 

o22   =   —  sin  <p  sin  ^ 

+  cos  <f>  cos  ^  cos  8 

at32   =    —  sin  8  cos  \f> 

Zi 

<tl3   =    sin  e  sin^s 

a23    =   sin  8  cos  <f> 

«33     —     COS  8 

With  the  coordinate  system  defined,  the  equations  of 
motion  may  be  expressed  in  terms  of  the  six  generalized 
coordinates  z  ,  R,  #"  ,  9,  V   and  <f>   by  the  use  of  Lagrangian 
mechanics.   The  only  simplifying  assumption  made  in  the 
derivation  of  the  equations  is  that  the  airplane  develops 
zero  thrust  in  the  spin.   This  is  frequently  the  case  in  an 


Reference  2,  page  158 
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actual  spin  situation  due  to  jet  engine  stalls  and/or  flame- 
outs  resulting  from  inlet  airflow  distortion.   However,  the 
restriction  imposed  by  this  assumption  can  be  easily  removed 
by  including  negative  dissipative  work  terms  on  the  right 
hand  sides  of  the  Lagrange  equations  which  follow.   Such 
terms  would  include  X,  Y  and  Z  force  components  as  well  as 
force  and  moment  derivatives  about  all  three  axes  due  to  o( 
and  {3   ,  o(.   and  3 ,    and  the  angular  rates  p,  q,  and  r. 

The  derivation  of  the  equations  of  motion  is  best  illus- 
trated by  first  writing  the  general  Lagrange  equation: 


(1) 


Where : 

T  =  System  kinetic  energy 

V  =  System  potential  energy 

q  =  A  generalized  coordinate,  r  =  l,2,---,6 

Fq   =  A  generalized  force  due  to  non-conservative 
forces,  r  -  1,2, •••,6 

For  this  particular  coordinate  system,  T  can  be  written  as: 

Where  I  I  I  ,  Cd x#  60 y/  &) z»  are  moments  of  inertia  and 
angular  rates  of  rotation  about  the  airplane's  principal 
axes. 

Using  Figs.  1  and  2,  the  angular  rates  relative  to  the 
inertial  reference  can  be  written  in  terms  of  the  spin 
coordinates  as  follows  : 
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fifc*  (£*■*)**  j  +ocos4>  (3) 


cjy  =  (<i>+i)o(zt  -esiup  <4> 


6^  =(1^+^)^3+47  (5) 

Note  that  <f  is  included  to  account  for  the  rotation  of  the 

X, ,    Y. ,    Z.    system. 

_  •         ••••• 

The  partial  derivatives  ^T/^  for  ^.r.:z^4j^j^jOJ^J 

and  0  are : 

|T*M«D  (6) 


** 


^J  «M^V  + I*^X^3  +^03/^3  +  12:^^3 


(7) 


(8) 


jE  r  IxCOxCos^    -lyCAySltip  (9) 


2l  *  Ix6)x^+Iy6)y^  +1*  Wgo^  (10) 


$-!««•  U1) 
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And  the  time  derivatives  of  these  are: 


4.  {&L  \ ,  i 


(14) 


5?©"  £tf4*ij +<*-*$)  +  lr{&i<<2a  +<U/*z3) 


(16) 


Next,  the  partial  derivatives  ^T^j  are  obtained: 
ll^O  (18) 


|£  -  M3r2fc  (19) 


|I»0  (20) 
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And  from  V  -  MgZ   the  values  of  °/(    are 


(24) 


(25) 


Next,  evaluating  the  quantities    CO^   ,    oC\  3   and  ^  'YJa,, 

where  k  ~  x,  y  or  z  and  i  ~  1,  2 ,  or  3  ,  and  r  =  4,  5,  or  6: 

* 


tiisCV  +  Yjetst  +  CV+Vtjs  +4>  (28) 


*>3  =  4>co$esjM<t>  +4>sin  q  cq$4>  (29) 


£2£  =&  COSG Co$ <f>  -<pStKl  6  S/AJ&  (30) 
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c<3s  =  -GSiue 


(31) 


^O  (32) 

*2>.«ie<»*  (33) 


.pi?  _  e<?s  e  cos  <£ 


(34) 


|f--**ew#  .      (35) 

jM*  -  S/M  ©  (36) 


^F        ^       TV  ~  7?"  ; 


(37) 


The  partial  derivatives  «5<--^-*  /  etc.,  are  also  easily  obtained, 
and  will  not  be  shown  here. 

Finally,  by  substituting  the  derivatives  defined  in 
equations  (12)  through  (25) ,  six  Lagrange  equations  of  motion 
can  be  written  in  terms  of  the  quantities  defined  in  equations 
(26)  through  (37)  and  the  partial  derivatives  of  the  angular 
velocities  are  as  follows: 

M'4q    +M<]   =  F2o  (38) 
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Mf>-Mrzfi~^  (39) 

+X£  C^€°<33  +  6fe*j*J  -  />  (40) 

Ix(faCd*<l>  -4d*SM4)  -2y(<bySM  <ptcOy<pcos<l>) 

•~(W^+ x^^  ^^l?;-*        (4D 

XJr$k*i  J  +•  &>x  *k)  +  1/  fay  o<ij  4  C0^  ^2 j  ) 

-  +1*  (  <&g  *3j  +  6J€  ^3)  =  /^  (42) 

j^^(xx^^,  ixcoy^;  =  F^  (43) 

As  discussed  in  Ref.  2,  the  concept  of  virtual  work  can 
be  employed  to  derive  the  non-conservative  forces  p     ,  Fa    , 
/v  /  Fq  /  f~ur  1    and1  P^    which  constitute  the  generalized  forces 
in  the  Lagrange  equations  (38)  through  (43).   In  this 
technique,  each  of  the  generalized  coordinates  (Z  ,  R,  Y   , 
0,  ty   or  p    )  are  increased,  one  at  a  time  by  a  small  positive 
increment  while  holding  all  other  coordinates  fixed,  and  the 
expressions  for  the  virtual  work,  terms  are  determined.   Since 
system  forces  and  moments  due  to  aerodynamic  effects  are 
applied  at  the  eg  and  expressed  in  orthogonal  X,  Y  coordinates, 
no  coupling  of  the  virtual  work  terms  occurs  and  they  can  be 
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readily  derived  by  inspection  using  pigs.  1  and  2.   The 
partial  derivative  of  the  sum  of  the  virtual  work  terms  is 
then  taken  with  respect  to  one  of  the  generalized  coordinates 
and  the  corresponding  generalized  force  is  thus  obtained. 
For  example  : 


Where  c^  ,  &y    ,  and  £^(  are  positive  incremental  angular  dis- 
placements about  the  respective  axes  X, ,  Y, ,  and  Z, . 
Since,  by  inspection  of   ig.  1: 

■££'=-/  (45) 

And,    as    a    result    of   orthogonality: 

^t.i/Lr  "5£x,  _  5|y,  ,  3§a   .  ^  (46) 

3»T         C7^*  <?£<>  J^  <*£* 

Thus  we  obtain: 

ft.=  -ft,  <47> 

And,  by  identical  procedures: 

F&=-Fxt  (48) 


*—* /*■*** 


(49) 


Fe-Mk  CoSip-tM^SttJ  </>  «  MKCdS4>+MyM<p  (5o) 
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+  Mi,  cos  e   =  M%  (51) 


fy   =  #2, 


(52) 


Now,  by  employing  the  direction  cosine  relationships  listed 
in  Table  I  these  forces  and  moments  can  be  obtained  in  terms 
of  forces  and  moments  along  and  about  the  airplane's  principal 
axes.   Thus  equations  (47)  through  (52)  can  be  rewritten  as: 


%o  =  *  (°<iiFk  +o(n/:y  +*S3F*) 


fif  -kif*  +  *llFY+*S/':*) 


(53) 


(54) 


F}    ~    -%fatlP*^2lFY  +  *3lFz)+(*/3ttxMllHy+°(llMl)       (55 


Fe-   Afx  COS<f>  -  MyS/Mtp  <56> 


Fp    -    M£  (57) 
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F*  =  (*,$  Mx  +°<ZS  fy  *  *33  Mi  )  <58) 


f 


Finally  the  complete  dynamic  equations  governing  motion  in  the 
chosen  coordinate  system  can  be  written  in  terms  of  the  six 
generalized  coordinates  z  ,  R,  if,  Q ,  ty ',   and  y       by  equating 
the  left  hand  terms  of  equations  (38)  through  (43)  with  the 
respective  right  hand  terms  of  equations  (53)  through  (58) . 
These  equations  are  written  in  their  entirety  in  Appendix  A. 
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IV.   INCORPORATION  OF  CONVENTIONAL  AERODYNAMIC  DATA 

A.   DERIVATION  OF  REFERENCE  PARAMETERS 

In  order  to  utilize  conventional  aerodynamic  force  and 
moment  coefficient  data  in  the  equations  of  motion  just 
derived  it  is  necessary  to  obtain  such  typical  reference 
parameters  as  airplane  angle  of  attack  (<K) ,    sideslip  angle 
{(3 )  /  center  of  gravity  velocity  (  V^a   )  *  r°H  rate  (  p  )  , 
pitch  rate  ( Q  )  ,    and  yaw  rate  (  f*  )  in  terms  of  equation 
variables  and  their  time  derivatives.   This  can  be  done  with 
the  aid  of  Figs.  1  and  3  as  follows: 

Vc9=  */?,  +y,  e^^Cs  (59) 

where  2.  ,  C?2  ,  and  Cj  are  unit  vectors  along  the  x,  /  Y-,  / 
Z ,  axes  respectively. 
But  from  Fig.  1: 

tf,  -  -A  ,  i  a  -to  i>  '-  -z° 

Thus  : 

l/c^-fce,  -ft?*-2;?j  (6o) 

Therefore  the  relative  wind  at  the  airplane  eg  is: 

Vaw~  *  ^  +tRSj,+itWs  (6i) 

By  use  of  the  direction  cosines  given  in  Table  I,  the  relative 
wind  can  be  resolved  into  components  along  the  airplane's 
X,  Y  and  Z  axes  as  follows: 

Va»x -  <*„  R  +  *,z  f*  +  °<>3  K  (62 ' 
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(63) 


(64) 


X 


« 


V 


RWZ 


*  Y 


tan<P^~ 


V 


RWZ 


V. 


V 
tan  {3  =   . 


RWY 


RWX  RWX 

Fig.  3   Definition  of  Positive  &(   and  jQ    Angles 
Therefore  o(    and  (3    can  be  written  as: 


oi^T&u     ( — — tz ; — ; — J 


(65) 


(66) 


The  magnitude  of  the  relative  wind  can  be  written  from 
equation  (61)  as: 


Vm*  {il*  CWz+z,zj 


V* 


(67) 
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Finally,  the  rotational  rates  p  ,  q  ,  and  r  about  the  air- 
plane's  x,  y  ,  and  z  axes  are  identical  to  the  angular 
rates  (jd     ,  COy   an<3  CO*   given  by  equations  (3j  ,  (4),  and 
(5),  and  can  be  written  directly  as: 

(68) 


P^( (p  +rJc//3  f-  e cos<p 


%*(p+f)^z$-e$tti<p  (69) 

r^  (<j>+V°<33  +  4  (70) 


B.   MODELING  THE  FORCES  FOR  THE  STEADY  SPIN 

Tabulated  aerodynamic  data,  although  virtually  always 
appearing  in  the  form  of  dimensionless  coefficients,  varies 
widely  in  terms  of  the  reference  axis  system,  point  of 
application  of  forces  and  moments,  and  the  characteristic 
length  used  in  non-dimensionalization,  depending  on  the 
intended  application  of  the  data.   By  the  use  of  the  princi- 
pal axis  system  in  the  foregoing  derivation  of  the  equations 
of  motion,  considerable  simplication  was  achieved.   There- 
fore it  appears  that  conversion  of  aerodynamic  data  to  the 
principal  axis  system  and  eg  as  references  would  generally 
be  preferable  to  the  alternative  of  adapting  the  equations 
to  the  particular  form  of  the  data.   This  should  present  no 
serious  problems,  however,  since  the  equations  must  be 
solved  numerically  and  conversion  of  the  data  by  computer 
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is  by  comparison  a  relatively  minor  additional  task.   Since 
some  data  manipulation  will  likely  be  required  in  any  event, 
the  mean  aerodynamic  chord  ( C )  has  been  used  in  the  follow- 
ing development  as  a  common  characteristic  length  in  an 
attempt  to  be  consistent  and  also  allow  the  use  of  a  single 
"characteristic"  length  parameter  as  a  computer  program 
input.   However,  the  equations  which  follow  could  be  adjusted 
appropriately  to  use  other  characteristic  lengths. 

The  general  form  of  the  aerodynamic  forces  and  monents 
contained  in  the  right  hand  terms  of  equations  (53)  through 
(58)  can  be  expressed  as  follows: 

Fi  '  '/%  cF;  ?v*s     i  ■  x,  Yj  <"-  i 
Ml  =  hC^^Sl 

Where  C    has  been  used  as  a  characteristic  length  in  non- 

dimensionalizing  all  moment  coefficients.  The  force  and 

moment  coefficients  (  Cpi  and  Cf*,'  )  are,  in  turn,  of  the 
following  general  form: 

CFi  =  CF!l  +  CFil  ■  Hp)  +  CFJ3  ■  f  fo)+  CFiA  -f(r) 
+  Cfis--  f(U)  +  Cpiffdr)  +CFi7  -fdO 

Where  : 

1.      All  Cpjj    and  Cft  ;j     (  j  =    1 ,  2 ,  •  •  -  7    and    i    =  x,y   or    z) 

are  aerodynamic  force  and  moment  derivatives,  as 
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defined  in  Table  II  below,  which  must  be  experi- 
mentally determined  and  tabulated  for  a  wide  range 
of  pi  and  A   values. 

Table  II :   Definitions  of  Aerodynamic  Derivatives 


SUBSCRIPT    j 

i~x,y , 
or    z 

1 

2 

3 

4 

5 

6 

7 

cm 

%(&) 

3»CW 
3*« 

c»ij 

<*ifat) 

2.  t(P),   t(^)  ,  and  f(r)  are  functions  of  roll,  pitch 
and  yaw  rate  respectively. 

3.  £     ,  ^  and  ^  are  defined  as: 
oq.  ~   Aileron  deflection 

&r   =  Rudder  deflection 
Sh  ~   Elevation  deflection 
A  good  example  of  the  form  of  tabular  data  that  could  be  used 
in  this  manner  is  contained  in  Ref.  3. 

At  this  point  some  comments  are  appropriate  to  insure 
that  the  above  method  of  handling  the  aerodynamic  forces 
and  monents  are  put  in  proper  perspective.   First  of  all  it 
should  be  noted  that  the  validity  of  any  calculations  of 
spin  parameters  obtained  by  this  approach  to  the  spin 


Aerodynamic  derivatives  with  respect  to  oK   and  fS   have  not 
been  included  since  6(  and p   are  zero  for  the  steady  spin, 
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problem  depends  entirely  on  the  modeling  of  the  aerodynamic 
forces  and  monents  since  the  equations  of  motion  as  derived 
in  Section  III  are  exact.   Secondly  it  should  be  recognized 
that  the  aerodynamics  of  the  spin  problem  involve  low  speed 
three  dimensional  stalled  flow,  and  as  such  forces  and  moments 
are  nonlinear.   Thus  the  superposition  of  the  separate 
effects  of  triaxial  rotation,  control  deflection,  and  other 
static  wind  tunnel  data  cannot  be  justified  analytically. 
The  only  wind  tunnel  test  technique  that  presently  appears 
to  offer  any  possibility  of  yielding  satisfactory  aerodynamic 
data  for  spin  is  the  rotary  balance  force  and  moment  measuring 
apparatus  which  is  reported  on  in  Ref.  4. 

Next,  allowing  that  for  want  of  better  data  the  conven- 
tional approach  outlined  above  will  be  used  to  exercise 
computer  program  and  look  for  a  zeroth  order  approximation 
of  spin  characteristics,  there  is  another  problem  that  stems 
from  this  application  of  conventional  data.   The  problem  is 
that  existing  data  is  generally  of  insufficient  range  in  o( 
and  /3  to  accommodate  the  high  angles  of  attack  and  side- 
slip angles  characteristic  of  autogyration.   Also  many  of 
the  basic  set  of  42  aerodynamic  coefficients  suggested  by 
the  above  development  describe  coupling  effects  which  in 
conventional  airplanes  at  controlled  flight  altitudes  are, 
in  general,  not  significant  and  therefore  have  not  been 
measured.   Thus,  for  example,  roll  moment  due  to  elevator 
deflection  is  normally  ignored,  but  may  be  of  some  signifi- 
cance in  a  spinning  F4  airplane  where  the  strong  negative 
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dihedral  of  the  stabalator  could  cause  the  upwind  panel  to 
act  as  a  rudder  with  prospin  deflection  at  high^  and/3 
angles.   Certainly  in  such  cases  none  of  the  aerodynamic 
derivatives  including  the  42  above  can  be  ignored  without 
testing. 

Now  admitting  to  the  fact  that  force  modeling  with  con- 
ventional stability  derivatives  will  be  done  with  at  best 
an  incomplete  set  of  data,  it  then  becomes  necessary  to 
exercise  caution  to  insure  that  those  coefficients  that  are 
available  and  used  will  allow  a  force  and  moment  balance  to 
be  somehow  achieved.   Thus,  if  at  least  one  damping  term  is 
not  included  for  every  driving  term  in  each  of  the  six 
degrees  of  freedom,  obviously  no  stabilized  solution  is 
attainable.   Finally,  in  fairness  to  this  effort  it  must  be 
pointed  out  that  the  modeling  of  aerodynamic  forces  and 
moments  on  the  spinning  airplane,  although  recognized  as 
the  key  to  the  problem  under  study,  was  non-the-less 
peripheral  to  the  expressed  purpose  of  this  preliminary 
work.   Once  a  workable  computer  program  has  been  written 
and  checked  by  use  of  conventional  aerodynamic  data,  sub- 
sequent efforts  should  be  directed  towards  improving  the 
modeling  of  aerodynamic  forces  and  moments. 
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V.   SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

A.   GENERAL  CONCEPT 

From  an  inspection  of  the  equations  of  motion  contained 
in  Appendix  A  (in  abbreviated  format)  it  is  evident  that, 
although  they  are  ordinary  second  order  differential 
equations,  they  are  also  so  extensively  coupled  and  non- 
linear  that  direct  analytical  integration  is  out  of  the 
question.   Even  if  the  left  hand  sides  of  the  equations 
were  linear,  the  forces  and  moments  which  constitute  the 
right  hand  sides  are  functions  of  a  three  dimensional  flow 
involving  the  geometry  of  the  airplane.   Consequently  any 
solution  scheme  will  of  necessity  involve  some  form  of 
iterative  numerical  procedure  to  circumvent  direct  inte- 
gration.  The  problem  thus  becomes  one  of  finding  an 
iterative  scheme  which  provides  some  degree  of  uniform  con- 
vergence to  steady  or  quasi-steady  values  of  the  six  spin 
parameters  (  £?^,  %.    ,    V  ,  0  ,  *f    ,    and  p  )  without  losing 
numerical  significance  or  using  excessive  computer  time. 

If  the  problem  is  visualized  as  a  trial  and  error 
positioning  and  orientation  problem  where  the  airplane  is 
incrementally  reoriented  angularly  in  response  to  an  aero- 
dynamic moment  inbalance  and  simultaneously  adjusted  in  ^  / 
/and  R.   in  response  to  an  aerodynamic  force  inbalance,  it 
then  appears  to  be  reducible  to  a  simple  series  of  stepwise 
adjustments  or  iterations.   Also,  since  the  problem  would 
thus  be  restarted  on  each  successive  iteration  with  updated 
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force  and  moment  data  corresponding  to  the  adjusted  spin 
parameters,  such  a  procedure  would  avoid  the  possibility  of 
loosing  numerical  significance.   The  problem  of  determining 
consistent  incremental  changes  in  each  spin  parameter  to 
achieve  some  degree  of  uniform  convergence  of  such  an 
iterative  scheme  can  be  easily  handled  by  using  the  principle 
of  conservation  of  momentum.   By  rearranging  the  equations 
of  motion  in  such  a  way  as  to  facilitate  the  computation  of 
accelerations  in  each  of  the  six  spin  coordinates  it  is  then 
possible  to  apply  a  common  increment  of  time  and  thereby 
compute  a  dynamically  consistent  set  of  incremental  changes 
in  the  parameters  and  obtain  new  parameter  values  for  the 
subsequent  iteration.   This,  in  a  sense,  amounts  to  follow- 
ing the  motion  of  the  airplane  as  is  done  in  the  conventional 
moving  coordinate  system  approach  to  the  problem,  but  the 
hybrid  inertial  coordinate  system  has  several  distinct 
advantages  in  this  application,  as  discussed  below. 

The  rationale  upon  which  this  suggested  solution  techni- 
que is  based  is  that  it  duplicates  in  a  computational  sense, 
what  actually  happens  to  an  airplane  under  such  flight  con- 
ditions.  However  one  key  advantage  that  the  computational 
scheme  has  over  the  actual  airplane  is  that,  on  each  restart 
or  iteration  we  can  discard  any  angular  momentum  inconsistent, 
with  a  steady  spin  situation.   Thus,  we  can  avoid  the  problem 
of  carrying  along  tumbling  momentum  which  could  considerably 
lengthen  the  computational  time  or  cause  the  solution  to 
pass  right  through  a  steady  spin  mode  or  recovery  without 
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stabilizing.   Of  course  to  realize  this  advantage  it  is 
necessary  to  be  able  to  differentiate  between  that  part 
of  the  motion  which  constitutes  tumbling  and  that  part  which 
is  consistent  with  steady  spinning.   Herein  lies  one  of  the 
major  advantages  of  writing  the  equations  of  motion  in  an 
inertial  cylindrical  reference  frame  rather  than  in  a  body- 
fixed  moving  coordinate  system,  for  in  the  hybrid  inertial 
system  tumbling  is  readily  recognizable  as  terms  involving 
e,ty      or0. 

B.   SIMPLIFICATION  OF  THE  EQUATIONS  OF  MOTION 

With  the  above  approach  in  mind  the  equations  of  motion 
can  now  be  simplified  sufficiently  to  meet  the  requirements 
of  the  solution  scheme  without  the  necessity  of  making  any 
restrictive  simplifying  assumptions.   Also  since  the 
intended  iterative  technique  deals  with  the  problem  directly 
in  dynamical  units,  we  do  not  have  to  non-dimensionalize  as 
would  be  necessary  if  one  resorted  to  a  somewhat  more 
abstract  mathematical  method  of  varying  the  spin  parameters 
to  obtain  convergence.   As  a  result  one  can  avoid  obscuring 
the  physical  nature  of  the  problem  and  is  in  a  better 
position  to  utilize  insight  in  monitoring,  trouble  shooting, 
and  interpreting  iteration  trends. 

The  simplification  is  started  by  recognizing  that  only 
steady  or  quasi-steady  spin  modes  are  sought.   Thus,  by 
setting  angular  orientation  rates  (  Q,    tp   ,  and  <p  )    equal  to 
zero  a  considerable  number  of  terms  can  at  once  be  eliminated. 
The  equations  contained  in  Appendix  A  thus  reduce  to  the 

following  form: 
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^Equation: 

m£6  +  m$  =  -  (cF*  «„  +  eF,  ofz3  +  cFi  »c33)  to  {(tt)  * 


\  Equation: 


MR-Mt'R.  =  -  (CPx  of„  +CFy  c<ll  +  CFe  of  J  %{(tg)* 

*(*)**■&)*}  (74) 


V  Equation: 

mr  (zni+Rif)+ix{(P+lr)t(,3  +  ecosip}^  -tiy{(f 

+?)*<23-  &  StH<t>]  <<23  1- 1 i  {(  WJ  M3  j  ■/■  $>}  «J3 


(75) 
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^Equation : 

-  esin<p}siti<p  -i*rzco$e  siN<t>dl3  (?g) 

-ZytZC0S  &  COS  <f>  o{23 -h  X^z  S/A/&  o(33 

\U  Equation : 

+lBJ($+*)'<>i+0fai  (77) 

0Equation : 

Although  R.    is  a  rate  parameter  which  is  not  consistent 
with  the  steady  spin  assumption,  it  effects  both  the  values 
of  c\  and  fo   and  also  the  magnitude  of  the  relative  wind ,  and 
therefore  must  be  retained  in  the  equations. 

Next,  to  put  the  equations  into  a  useable  form  for  the 
suggested  iteration  scheme  we  simply  solve  each  equation 
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for  the  corresponding  acceleration  parameter  in  terms  of  the 
other  variables.   The  equations  then  become: 
Z    Equation: 

K  Equation: 

M~tZR-  (W,,  +  CFf  th+G^&U  W  (80 ) 

}f  Equation : 

-(Xx-ly)eStU6  Slti<t>  C0S<f>  -Z{Qc<n-2M&ti  (81) 

Q  Equation : 

0  -  [l,  COSl(p  +  ly  S/M  Z<P  }    [-(l%-Iy)($+  H)  SM  Q  Slti  j>  COS  (f) 
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I//  Equation: 

-I 


(83) 


~(lx-lr)OMe  Siti<pCOS<t)-Iz  0COS  &+(CM*e<ts 


(j)  Equation: 

$=(•&£&) y*-M9stfi><b cos <p  ~((f+r)cose  (84) 


Before  proceeding  to  develop  the  details  of  the  proposed 
iterative  scheme  of  solving  these  equations,  some  comments  on 
the  steady  spin  assumption  should  be  made.   First  of  all, 
truly  oscillatory  spin  modes,  as  opposed  to  lightly  damped 
steady  spin  modes  that  may  be  very  slow  in  convergence 
should  not,  in  general,  result  from  the  use  of  mean  force 
and  moment  model  test  date,  the  usual  form  of  such  data. 
Thus  if  the  oscillatory  nature  of  a  spin  were  attributable 
to  periodic  vortex  shedding  from  the  fuselage  or  stalled 
aerodynamic  surfaces  (von  Karman  Vortex  Street  phenomenon) , 
the  use  of  mean  aerodynamic  forces  and  moments  should 
effectively  filter  out  such  oscillation  from  the  computer 
solution.   If  however  the  oscillations  were  strictly  a 
result  of  the  dynamics  of  the  problem,  such  as  a  limit 
cycle  type  of  oscillation,  then  the  oscillatory  nature  of 
the  spin  should  appear  in  the  computations  as  a  failure  to 
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converge  to  steady  spin  parameter  values.   In  this  case,  the 
complete  equations  would  have  to  be  used  along  with  a  very 
small  time  increment  and  the  stability  derivatives  due  to  o( 
and  (3  to  obtain  an  accurate  description  of  the  motion  in  the 
form  of  a  numerical  time  history.   Although  this  would  involve 
solving  the  full  equations  for  the  corresponding  acceleration 
terms  as  was  done  in  equations  (79)  through  (84)  and  would 
also  lengthen  the  computation  time  it  would  not  be  a  diffi- 
cult extension  and  might  well  provide  some  valuable  insight 
into  the  stability  of  such  oscillatory  modes.   This  form  of 
the  spin  equations  would  be  especially  valuable  in  studying 
the  effect  of  applying  oscillatory  recovery  control  inputs  in 
resonance  with  the  spin  oscillations  to  obtain  an  "autopilot 
spin  recovery"  from  an  otherwise  irrecoverable  spin  mode. 
Another  use  of  the  full  equations  would  be  to  study  the 
stability  of  spin  solutions  obtained  by  other  numerical 
methods.   However  terms  involving  Q  ,  [b    ,    and  0  would  re- 
introduce the  effects  of  angular  momentum  into  the  computations 
and  would  probably  cause  excessive  tumbling  if  the  full 
equations  were  used  to  initially  solve  for  spin  modes.   There- 
fore either  the  following  computational  scheme  or  some  other 
numerical  technique  for  solving  non-linear  differential 
equations  such  as  the  iterative  procedure  described  on  page 
270  of  Ref.  7  should  be  employed  to  initially  solve  for  spin 
modes . 

C.   SUGGESTED  COMPUTER  SOLUTION  SCHEME 

Although  the  equations  are  rather  complex  the  suggested 

computer  iteration  scheme  is  fairly  simple.   Basically  some 
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starting  position  or  initial  "guess"  must  be  supplied  in 
terms  of  altitude  rate,  spin  radius  ,  spin  rate,  orientation 
(9,  jC  ,  and  <p  )  ,  along  with  control  positions  X  ,  C      and  (> 
Angle  of  attack  and  angle  of  sideslip  are  then  computed  from 
equations  (65)  and  (66)  and,  in  turn  C\   and  G   are  used  to 
obtain  values  for  the  various  aerodynamic  derivatives  dis- 
cussed in  Section  IV  B  by  interpolating  from  appropriate 
tables.   The  six  aerodynamic  force  and  moment  coefficients, 

Cf:    andC^f'f    are  then  computed  as  outlined  in  equations  (71) 
and  (72)  and,  with  these  values  in  hand,  ifd  and  A  can  be 
computed  directly  from  equations  (79)  and  (80) .   The  other 
accelerations  }f   ,  C?  ,  ft   ,  and  0   ,  can  be  computed  by  a 
looping  routine  encompassing  equations  (81)  through  (84) . 
This  looping  will  be  necessary  due  to  the  coupling  of  the 
respective  acceleration  terms  in  each  of  these  equations. 
Once  consistent  values  of  y  ,  Q   ,  lp    ,  and  <£)    have  been 
obtained  from  the  looping  iteration,  incremental  changes  in 
each  of  the  seven  computation  parameters  can  be  obtained 
from  the  following  equations: 

£,/   >*i   +  £<3t  (85) 

Am***^*4*  (86> 


^r%+k^fulT  <87> 


Radius  rate,  fc   ,  is  initially  set  equal  to  zero, 
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jU-V*  (88) 


%.Q    -  V®  *■  (89) 


^♦,;=«VHr  <90) 


hwrKi+b^T  (91) 

If  the  full  equations  were  to  be  used,  obviously  it 
would  then  be  necessary  to  also  compute  Z?  ,  If    ,  ty    ,  &  , 
(// ,    and  <p      as  follows: 


vv^ 


(92) 


£   v=^  -+&,       Ai-h&^i  (93) 


«%r*0  s  Pw  *  ^At 


(94) 


SVj^  +  W*^ 


',',  At* 


(95) 


^NsAw^4t  (96) 
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^VW^*^  (97) 


The  computational  sequence  outlined  above  is  then  repeated 
using  the  revised  values  of  the  parameters  involved  and 
continued  until  some  pre-designated  convergence  criteria  is 
met,  such  as  all  accelerations  being  less  than  some  small 
epsilon  value  or  the  rates  of  change  of  the  various 
parameters  are  each  in  turn  less  than  some  small  character- 
istic values  as  individually  specified.   For  example,  one 
foot  per  second  in  £  ,  0.5  feet  per  second  in  £  ,  one  degree 
per  second  in  y  ,  etc. 

The  choice  of  an  incremental  time  probably  will  have  a 
very  significant  effect  on  the  overall  success  of  this 
iteration  scheme.   The  use  of  a  very  small  value  of  ^'twill 
of  course  provide  the  smoothest  convergence  but  at  a 
significant  expense  in  computer  time.   A  large  £t>    on  the 
other  hand  can  easily  preclude  convergence  by  drastically 
altering  all  parameters  from  iteration  to  iteration.   To 
illustrate  this  point  and  also  pin  down  what  is  meant  by 
large  and  small  A*C  it  is  necessary  to  step  back  a  bit  and 
consider  conceptually  what  is  really  being  done  by  using  an 
iterative  solution  scheme.   As  pointed  out  earlier,  the 
equations  are  highly  non-linear  and  extensively  coupled, 
and  consequently  one  has  very  little  insight  as  to  what 
linearizing  approximations  would  be  valid  or,  for  that 
matter,  over  what  range  of  values  can  the  effects  of  an 
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individual  parameter  change  be  considered  linear.   This 
problem  is  however  overcome  in  the  suggested  iteration 
scheme  by  allowing  the  dynamics  of  the  problem  to  pick  the 
magnitude  of  the  incremental  changes  in  each  parameter. 
Thus  if  the  ^l»  used  is  of  such  a  magnitude  as  to  constitute 
a  small  pertubation  relative  to  the  characteristic  time  of 
the  problem,  the  time  per  revolution,  we  essentially  estab- 
lish appropriately  sized  small  pertubations  in  each  of  the 
spin  parameters  as  well.   Thus,  for  example,  a  spin  rate  of 
40  per  second  would  yield  a  9  second  turn  and  possibly  a 
good  starting  guess  for  an  appropriately  sized  ^"t  would  be 
.09  seconds.   The  validity  of  this  concept  will  not  be 
rigorously  or  otherwise  established  here,  other  than  to  say 
that  in  general  small  pertubation  techniques  have  enjoyed 
some  success  in  other  engineering  applications  and  this 
appears  to  be  a  quite  natural  extension. 

The  problem  of  picking  appropriate  starting  guesses  for 

2L  ,  r\   t   ¥    >   0  ,   Y    >    anc^  *P    unfortunately  is  not  as  easily 
handled  as  deciding  on  an  initialZrt.   Although  some  rather 
broad  guidelines  can  be  stated  such  as  steep  nose  down 
pitch,  large  radius  and  low  spin  rate  for  normal  erect 

(upright)  spins  and  shallow  nose  down  pitch,  small  radius 
and  high  spin  rate  for  flat  spins,  the  actual  searching  for 
individual  modes  must  proceed  on  a  trial  and  error  basis 
until  an  analysis  of  solutions  suggests  some  better  pro- 
cedures.  Starting  guesses  that  don't  lead  to  spins  will  of 
course  show  up  as  recoveries  to  linear  stalled  flight  atti- 
tudes or  complete  fly-away  recoveries. 
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Certainly  some  insight  can  be  used  to  pick  Z^  as  a 
negative  (down)  and  of  an  order  of  magnitude  consistent 
with  stalled  flight  velocities  (T.A.S.)  at  the  altitude 
where  the  solution  is  sought.   Likewise,  a  wings  level 
attitude  with  the  nose  down  and  pointed  at  the  spin  axis, 
e.g.,  9  =  +30°,  P  =   -90°  and  P  =  +90°,  would  be  a  reason- 
able starting  orientation  for  an  erect  right  spin.   Finally, 
%   should  obviously  be  chosen  small  so  that  it  will  not  toss 
the  airplane  out  of  the  ball  park  under  the  action  of 
centrifugal  force  before  the  aerodynamics  can  come  into 
play.   One  point  should  be  made  however  concerning  searches 
for  flat  spin  modes  with  small  radius  values.   If  conven- 
tional force  modeling  is  used,  as  discussed  in  Section  IV  B, 
then  in  addition  to  the  general  lack  of  validity  of  adding 
non-linear  aerodynamic  forces  and  moments,  there  is  the 
compounding  of  errors  due  to  the  fact  that  free  stream  o(  , 

y     and  &    vary  from  point  to  point  over  the  fuselage  and 
wings  as  a  function  of  orientation  when  radius  is  reduced  to 
the  same  order  of  magnitude  as  fuselage  length  or  wing  span. 
Thus,  for  example,  if  the  wings  are  level  and  the  nose  is 
pointed  at  the  spin  axes  in  a  fairly  flat  attitude  the 
geometry  of  the  problem  will  yield  a  significant  difference 
in  both  sideslip  velocity,  and  consequently  /3  ,    from  nose  to 
tail  for  a  small  spin  radius.   Of  course  aerodynamic  data 
for  this  type  of  flow  condition  cannot  be  obtained  in  a 
static  wind  tunnel  setup. 
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D.   AUXILIARY  ANGULAR  RELATIONSHIPS 

The  particular  set  of  "Euler  Angles"  used  in  the 
derivation  of  equations  (79)  through  (84)  is  just  one  of 
several  variations  of  such  positioning  angles  which  are 
generally  classified  as  Euler  Angles.   The  principal 
advantage  gained  by  the  use  of  such  sets  of  positioning 
angles  is  that  they  describe  a  single  fixed  spatial 
orientation  of  a  rigid  body  for  given  values  of  Q>  ,  f    , 
and  (p   ,  and  yield  a  convenient  set  of  equations  for  the 
direction  cosines  with  the  order  of  angular  rotation 
implicit  in  the  equations.   However  unlike  the  more  con- 
ventional pitch,  roll  and  yaw  orientation  angles  which  can 
be  readily  visualized  once  the  order  is  specified , 
orientation  in  terms  of  Euler  Angles  is  difficult  to  visualize 
without  some  sketching,  which  can  be  readily  appreciated  by 
refering  to   Fig.   2.   Thus  unless  the  orientation  is  close 
to  quadrantal  a  computer  output  of  orientation  in  Euler 
Angles  may  be  difficult  to  interpret. 

To  alleviate  this  situation,  relationships  between  an 
ordered  set  of  orientation  angles  (1)  pitch,  (2)  roll,  and 
(3)  yaw,  and  the  Euler  angles  shown  in   Fig.   2  were 
derived.   The  derivation  was  facilitated  by  the  sketch  of 
the  spherical  triangle  relationships  between  the  two 
angular  sets  shown  in  Fig.   4  and  involves  the  use  of 
Napier's  analogies  from  Spherical  Trigonometry,  Ref.  5. 

The  use  of  Napier's  analogies  avoids  the  quadrantal 
ambiguity  associated  with  the  solution  of  right  spherical 
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triangles  by  other  methods.   Thus  the  equations  for  con- 
verting the  ordered  angular  rotations  Q> ,  $     and  </>     to  the 
set  of  Euler  angles  used  in  this  paper  are: 

-T..-rf*«S2l  1 


*.  *  -w'/Wfyr.,  #;j  -t.  j  { s«g  J 


^«-,{s^,J-^(*,®|5"*;} 


u-'faM[w'f'-M(%)Tw(g')}] 


0-27* 


(98) 


(99) 


(100) 


Fig.  4   Schematic  of  Auxiliary  Angular  Relationships 
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The  equations  for  obtaining  Qp  ,   0^     and  Yy      f rom  Euler  angl 
are  : 


es 


(101) 


.   a  .wfat&wCa/  ..-MP^Cfefr . 


fw  (t2±^ 


c^5  (*££«j 


>  :io2) 


f*-f<hZT<M 


(103) 


Equations  (101)  through  (103)  can  be  converted  for  computa- 
tions in  radians  by  substituting  'j^£  for  90°. 
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VI.   THE  INVERSE  PROBLEM 

As  discussed  in  Section  III  B  the  conventional  method  of 
modeling  the  aerodynamic  forces  and  moments  on  the  spinning 
airplane  with  static  wind  tunnel  data  lacks  mathematical 
validity.   Therefore  it  will  be  necessary  to  develop  a 
better  method  of  modeling  the  aerodynamics  of  the  problem 
if  the  computational  approach  to  predicting  spin  modes  is  to 
be  successful.   The  work  done  with  the  rotary  balance,  Ref. 
4  has  been  one  effort  in  this  area.   Such  work,  however,  is 
difficult  to  pursue  without  some  method  of  obtaining  the 
actual  forces  and  moments  as  a  basis  for  comparison. 
Apparently  in  the  past  the  validity  of  such  aerodynamic 
data  has  been  checked  by  the  indirect  method  of  using  it  to 
solve  the  partially  linearized  equation  and  then  comparing 
computed  spin  modes  with  known  spin  characteristics.   This 
approach  may  allow  some  estimate  of  the  overall  success  of 
the  combined  result  of  the  equations  and  data  but  doesn't 
establish  the  validity  of  either.   A  better  approach  would 
be  to  use  the  exact  equations  derived  in  Section  V  in  such 
a  manner  as  to  allow  comparison  of  the  6  computed  aero- 
dynamic forces  and  moments  with  the  corresponding  forces 
and  moments  on  the  actual  spinning  airplane.   Fortunately 
this  can  be  done  for  the  steady  spin  case  and  thus  allow 
some  progress  to  be  made  in  modeling  steady  or  mean  values 
of  forces  and  moments. 
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In  order  to  use  the  equations  of  motion  derived  in 
Section  V  in  this  application  it  is  first  necessary  to 
obtain  the  values  of  £tf ,  K    ,  &  ,  ft    ,  and  y    for  the  actual 
airplane  or  free  spinning  model.   For  the  case  of  the  model 
in  the  free-spinning  wind  tunnel  discussed  in  Ref.  6, 
determination  of  the  above  spin  parameters  can  obviously  be 
obtained  from  tunnel  speed,  and  photographic  coverage,  and 
thus  will  not  be  discussed  further.   For  the  actual  airplane 
the  means  of  getting  these  parameters  is  still  fairly 
simple  but  in  some  cases  less  obvious,  and  consequently 
will  be  outlined  here.   First  of  all,  Z-j   can  be  determined 
by  the  pilot  altitude/time  observation  or  photo  panel  data, 
or  externally  by  ground  station  precision  radar  or  optical 
tracking.   Likewise  J    can  be  determined  either  by  the  pilot 
or  from  ground  station  photographic  coverage.   Spin  radius 
can  be  obtained  by  resolving  triaxial  eg  accelerometer  data 
along  the  principal  axes  and  using  the  following  relationships 

■3?*3;*ti-*°+(t*tf  <104) 

Airplane  orientation  in  terms  of  Euler  Angles  &  ,   y    ,    and 
(D  can  be  determined  from  values  of  pitch,  roll,  and  yaw 
rates  and  the  triaxial  accelerations  as  follows : 


P*  isiue  $iN(p 


(106) 
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q  =  /SMS  COS<p  (107) 

r^fcose  (108) 

Thus  : 

([>=T<lHl{ 1 j  (109) 


S/N0 


and  from  Table  I : 
or  solving  for  ^  : 


(no) 


(111) 


Obviously  there  will  be  some  loss  of  accuracy  in  the 
determination  of  spin  radius  and  orientation  if  the  spin  is 
not  exactly  steady.   However  small  fluctuations  in  the 
triaxial  accelerations  and  pitch  roll  and  yaw  rates  can  be 
averaged  out  and  should  allow  computation  of  representative 
values  of  /^  ,  Q   ,    if/  ,    and  (h    for  the  near  steady  spin  case. 

Next,  equations  (79)  through  (84)  can  be  adapted  to  the 
problem  at  hand  by  setting  R  and  all  acceleration  parameters 
equal  to  zero.   The  equations  thus  reduce  to  the  following: 


50 


Z.   Equation: 

"  Mm^J} =  Cf*  *13  * cpy  «xs + ^  "*       (112) 

£  Equation : 

=     ^5f</+4>**/"*Q%*3/  (113) 


zy2rm 


0    Equation : 


Z(Xx-ly)YZS/AJZe  SM<t>COS<fi     _  ^ 


(114) 


0  Equation: 
ty  Equation : 

f%^^^a*V»S°  U16) 

<P  Equation : 


(117) 


Equations  (112)  through  (117)  can  also  be  cast  into  matrix 
format  for  direct  computer  computation  as  follows: 
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OXS0  ~S7/^     O 
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(119) 


^  (  Ix smz<p+1y  Cosz4>  -I?)fzsiue cose 
esc  {(H?+(.i.f} 
o 

esc{(t/i)z+(i6)3-} 


where  the  fact  that  Cu-s.O     ,    as  expressed  by  equation  (116), 

has  been  used  in  equation  (114)  to  show  that  Cpw=Q    in  the 

Yi 
steady  spin. 

In  addition  to  allowing  the  computation  of  the  force  and 

moment  coefficients  from  steady  spin  parameters  the  results 

of  this  development  show  that  Cy    ,  C/^v    and  C/^   are  not 
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necessarily  zero  in  the  steady  spin  case,  contrary  to  what 
one  might  have  intuitively  assumed.   In  fact,  the  steady 
spin  situation  where  C^    >C^     ,    an<3Lyu_  are  all  zero  is  a 
special  case  not  likely  to  be  encountered  in  airplanes  of 
conventional  design.   To  illustrate  this  point,  equations 
(115)  and  (117)  can  be  solved  for  the  case  where  the  moment 
coefficients  are  all  zero  with  the  result  that: 

l*SltiZ<p+  If  COsty^I^  (120) 

and 

&  -ly  (121) 

or  simply: 

Thus  it  is  seen  that  in  a  conventional  airplane  where  it  is 
unrealistic  to  have  equation  (122)  satisfied,  aerodynamic 
moments  must  be  present  to  balance  the  gyroscopic  moments 
and  gyroscopically  precess  the  inertial  orientation  of  the 
airplane  sufficiently  to  maintain  a  fixed  orientation 
relative  to  the  spin  axis  and  radius  vector. 
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VII.   RESULTS  AND  DISCUSSION 

A.   ANALYTICAL  ASPECTS 

Although  the  primary  intent  of  this  effort  was  to 

develop  another  method  of  computing  airplane  spin  modes 

from  wind  tunnel  test  data,  the  equations  of  motion  derived 

in  Section  IV  can  also  be  used  to  advantage  in  a  qualitative 

analytical  approach  to  the  problem  of  determining  spin 

characteristics.   In  addition  to  possibly  suggesting  new 

parametric  indices  of  spin  characteristics,  these  equations 

should  also  allow  a  number  of  experimentally  observed  effects 

of  parameter  variations  to  be  put  on  a  firmer  analytical 

foundation,  such  as  the  effect  of  the  sign  of  (I  -I  )  on 

3       x   y 

spin  recovery  [Ref .  8] .   The  fact  that  the  inertia   character- 
istics appear  so  prominently  in  equations  (79)  through  (84) 
confirms  test  results  which  have  suggested  that  the  gyro- 
scopic aspects  of  the  problem  are  as  important  as  the 
aerodynamics.   Another  result  of  this  analysis  of  the  spin 
problem  was  to  emphasize  the  non-linearity  of  the  aero- 
dynamics associated  with  the  spin.   Thus  computational 
techniques  which  use  superposition  in  modeling  the  aero- 
dynamics or  employ  partially  linearized  equations  of  motion 
should  be  expected  to  yield  rather  poor  results  where  high 
spin  rates  or  small  spin  radii  are  involved.   These  are 
probably  the  main  reasons  why  the  conventional  approach  has 
not  really  been  successful  in  treating  spins  in  general, 
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and  has  been  particularly  unproductive  in  predicting  flat 
spin  modes. 

B.   THE  COMPUTER  PROGRAM 

Although  the  suggested  computer  program  outlined  in 
Section  V  C  uses  conventional  aerodynamic  modeling  which  is 
not  entirely  representative  in  the  spin,  it  is  felt  that 
development  of  a  computer  program  capable  of  solving  the 
equations  of  motion  is  the  next  essential  step  in  this 
approach  to  the  spin  problem.   Once  the  computer  program 
is  working,  improvements  in  the  modeling  of  the  forces  and 
mements  can  be  incorporated  as  they  are  devised.   Thus  the 
computational  and  aerodynamic  parts  of  the  problem  can  be 
pursued,  in  the  interim,  as  separate  problems.   It  is  with 
this  in  mind  that  the  author  has  chosen  to  discuss  a  few  of 
the  aspects  of  the  computer  solution  problem  rather  than 
address  the  aerodynamics  of  the  problem. 

The  iterative  scheme  outlined  in  Section  V  B  may  exhibit 
a  tendency  toward  a  lightly  damped  oscillatory  form  of  con- 
vergence since  it  is  essentially  following  the  dynamical 
motion  of  the  airplane,  which  is  characteristically  lightly 
damped  in  a  spin.   The  oscillatory  nature  of  the  convergence 
is  compounded  somewhat  because  we  cannot  be  sure  how  far  a 
starting  guess  may  be  from  a  force  and  moment  equilibrium. 
To  improve  this  situation  there  are  a  number  of  ways  to 
provide  some  synthetic  damping  to  the  iteration  to  speed 
convergence.   Some  of  these  ways  are: 
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a.   If  excursions  in  radius  prove  to  be  a  problem  as 
a  result  of  momentum  in  radius  rate  (^) ,  the  motion  can 
be  damped  as  follows; 


'/ 


*o*/)=£*<»;4*At  (123) 

where  *3  is  some  factor  between  zero  and  one. 

b.   If  excursions  in  altitude  rate  become  a  problem  they 
can  be  handled  as  follows; 


*(n»)~*to)*$iAt  <127) 

Another  problem  that  may  arise  is  that  convergence  in 
orientation  may  be  deadbeat,  depending  on  the  size  of  <A"t  . 
This  is  due  to  the  fact  that  we  have  eliminated  orientation 
rates  to  reduce  any  tumbling  tendencies.   Thus  it  may  be 
necessary  to  augment  the  rate  of  change  in  orientation  per 
iteration  to  speed  convergence.   This  can  be  done  in  the 
following  manner: 


<w  <•*»+*(  M£  J 


(125) 


"    -L* 


Where  b  is  some  factor  greater  than  one. 

Finally,  since  our  starting  guess  in  Y  may  be  out  of  the 
ball  park  and  cause  large  changes  in  radius  prior  to  reaching 
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a  magnitude  consistent  with  the  aerodynamics  of  the  problem, 
it  may  be  necessary  to  augment  the  rate  of  change  of  this 
characteristically  unresponsive  parameter.   This  may  be 
accomplished  in  the  following  way: 


*(«>,)  -  &  +  c  (  ^ 


(128) 


where  C  is  some  factor  greater  than  one. 

C.   FOLLOW-ON  STUDIES 

As  mentioned  earlier,  the  equations  of  motion  as  derived 
in  this  paper  should  make  it  possible  to  provide  an  analyti- 
cal basis  for  many  of  the  heretofore  experimentally  observed 
characteristics  of  airplane  spinning  motion.   Such  a  cor- 
relation study  would  be  relatively  inexpensive  to  pursue  and 
would  probably  provide  a  wealth  of  additional  insight  into 
the  spin  problem  as  a  whole  and  possibly  be  beneficial  to 
the  more  immediate  goals  of  computer  program  development  and 
force-moment  modeling.   It  might  also  be  instructive  to 
program  the  full  equations  of  motion,  as  contained  in  the 
appendix,  on  an  analog  computer  having  a  digital  computer 
interface  capability  for  the  aerodynamic  data.   Some  other 
aspects  of  the  spin  problem  which  could  be  pursued  once  the 
computer  programming  and  force  modeling  problems  are  solved 
are  : 

(1)  Determine  the  spin  modes  of  currently  operated 
tactical  military  airplanes  and  investigate  the  effective- 
ness of  recommended  spin  recovery  techniques. 
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(2)  Determine  what  characteristic  aerodynamic  parameters 
can  be  most  economically  fixed  in  existing  airplanes  with 
irrecoverable  spins  and  investigate  the  magnitudes  of  such 
changes  that  would  be  required  to  make  recovery  possible. 

An  example  would  be  providing  rudder  stop  over-rides  or 
additional  deflection  for  spin  recoveries. 

(3)  Study  the  feasibility  of  providing  strap-on  devices 
such  as  a  deployable  fabric  ventral  fins  or  wing  mounted 
retro-rocket  pods  for  use  on  tactics  training  flights.   The 
size  of  the  device  or  necessary  thrust  requirements  could 

be  fairly  easily  determined  with  the  assistance  of  a  working 
computer  program,  or,  in  the  case  of  the  retro-rocket,  could 
also  be  determined  by  the  "Inverse  Problem"  Technique  out- 
lined in  Section  VI. 

(4)  Determine  "design  envelopes"  that  prescribe  limits 
in  terms  of  common  aerodynamic  derivatives  in  order  to 
assist  the  designer  in  avoiding  troublesome  combinations 
wherever  feasible.   Such  work  would  also  provide  the  basis 
for  improvements  in  the  current  Military  Specification 
MIL  F-8785  (ASG) . 

(5)  Study  the  feasibility  of  utilizing  an  autopilot 
to  apply  oscillatory  control  inputs  in  resonance  with  the 
dynamics  of  the  spin  to  break  an  otherwise  irrecoverable 
mode. 

These  are  only  a  few  applications  that  might  be  under- 
taken to  alleviate  the  spin  problems  inherent  in  the 
designs  of  todays  high  performance  tactical  military 
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airplanes.   However  the  satisfactory  accomplishment  of  just 
a  few  is  all  that  is  really  needed. 
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VIII.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  approach  to  the  spin  problem  outlined  in  this  paper 
should  provide  a  means  of  obtaining  some  additional  insight 
into  the  problem  of  aircraft  spin.   Once  a  successful 
computer  program  has  been  written  and  a  more  representative 
means  of  modeling  aerodynamic  forces  and  moments  has  been 
devised  it  should  be  possible  to  solve  the  equations  of 
motion  and  obtain  some  idea  of  the  spin  modes  that  a  given 
airplane  design  may  exhibit.   The  results  of  Section  VI 
indicate  that  gyroscopic  moments  play  an  important  role  in 
the  spin  problem.   Thus  it  may  also  be  worthwhile  to  treat 
the  spinning  airplane  as  a  slow  turning  free  gyroscope,  of 
relatively  high  inertia,  being  acted  upon  by  triaxial  aero- 
dynamic forces  and  moments.   Such  a  formulation  of  the  spin 
would  allow  the  application  of  the  broad  scope  of  knowl- 
edge of  gyroscopic  motion  from  the  fields  of  physics  and 
mechanics  to  the  airplane  spin  problem.   In  this  manner  it 
may  be  possible  to  obtain  some  additional  understanding  of 
the  stability  of  various  spin  modes  and  possibly  avoid  some 
undesirable  airplane  spin  characteristics  without  having 
detailed  knowledge  of  the  aerodynamics  involved. 
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APPENDIX 

SPIN  EQUATIONS 
The  following  equations  constitute  the  full  equations  of 
motion  for  a  spinning  airplane  as  written  in  a  cylindrical 
coordinate  system. 
£   Equation 


/v  Equation 
o    Equation 

+4>$\tiQ  cos<P) t  ecosfl - e <psi/o<f>}°<i3-tlx  {(^V,3 

+  QC0S<f>}{®  COS  6  $»</>  +  j>  S/AJ  &  COS  <f  }  +  Iy{($+t)o(z3 

+((p+f)(ecose  cos(f>  -<psifj&  si^<p)  --6s//0<p 

-&4>CoS<p>)^(2l  +ly  (((p+f)^  -  6  f/'J  <£>}  (&COSG  COSp 

-  4>s,n  &  s,tJ  4> }  +  h  { ( $■* fUn  *  (W)  ("  *  stu  &) 
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Equation 


2x  {($ + ?)*/  3  +■  (</>+$)(&  cos  e  sm4>  +4>  s/fJ  ecas<t>) 

+  &COS0  -Bp  SlK)<t) }  COS0- lx  <p  {((p+i) <Jl3-t€> C0S<pj Slti <f> 

-iy{(V>+f)*ki  +  (p+?)(ecos6  COS<p-0S/tie  sttifi) 

-6  SlfJ<t>   -  6<pG0S4>}S'liJ  0-1/0  { (fat) oil 3 

'€>ritj0}cos4>  -Ix  ((!/>■+  ?)o(t3  +  a  cosfi)  {(p 
+f) coss  cos<p}  i  it  {(fi+lrjto  +?}{(<P+f)s,'J6} 

T     Equation 

•4  \(<P+?)o(n  +  (<P+f)(6 cose  w<p + i)sme cos<p)-t S coscp 

-€><pS/*}(f>  /«/3  +  I*  {(<i>+lf)°(i3  f6<jOi^][6C6ie>S/^<p 
+  0  SIN  6  COS0  }+'*y  {(0+^)^z3+(^)(6  W6  COS0 
-<j)S/rJ  6  SlN<{))  -6  S/N0-  6  0  C0S(f)}p<z3  +  h  {(P 
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0  Equati 


on 


+ecos0}f(<P+t)  jVaJ  6  cos  4)  -  e  s»i  0}  -  Jy  {(tjufj^ 
- 6s//o 4>\{- (f*Y) s/a) e  s/tf 0  - e cos <p } 


-(Cm*  )(4r  {(**)*  +  <*)*+&)*} 


63 


BIBLIOGRAPHY 


1.  Johnson,  A.  E.   A  Literature  Survey  of  the  Problem  of 
Aircraft  Spins,  Masters  Thesis,  United  States  Naval 
Postgraduate  School,  Monterey,  California,  September 
1971. 

2.  Wells,  D.  A.   "Theory  and  Problems  of  Lagrangian  Dynamics", 
Schaum's  Outline  Series,  McGraw-Hill  Book  Company,  1967. 

3.  National  Aeronautics  and  Space  Administration  TN  D-6117, 
Utilization  of  a  Fixed  Base  Simulator  to  Study  the  Stall 
and  Spin  Characteristics  of  Fighter  Airplanes,  by  F.  L. 
Moore,  E.  L.  Anglin,  M.  S.  Adams,  P.  L.  Deal,  and  L.  H. 
Person,  Jr. ,  March  1971. 

4.  National  Advisory  Committee  for  Aeronautics  TN  2181,  The 
Aerodynamic  Forces  and  Moments  on  a  1/10  Scale  Model  of 
a  Fighter  Airplane  in  Spinning  Attitudes  as  Measured  on 
a  Rotary  Balance  in  the  Langley  20-Foot  Free  Spinning 
Tunnel,  by  R.  W.  Stone,  Jr.,  M.  Burk,  Jr.,  and  W.  Bihrle , 
Jr. ,  1950. 

5.  Kells,  L.  M. ,  Kern,  W.  F. ,  and  Bland,  J.  R.   Spherical 
Trigonometry,  McGraw-Hill  Book  Company,  1953. 

6.  National  Aeronautics  and  Space  Administration  Report 
R-57,  Status  of  Spin  Research  for  Recent  Airplane 
Designs ,  by  A.  I.  Neihouse,  W.  J.  Klinar,  and  S.  H. 
Scher,  1960. 

7.  Bekey,  G.  A.,  Karplus ,  W.  J.   Hybrid  Computation,  John 
Wiley  and  Sons,  Inc.,  1968. 

8.  National  Aeronautics  and  Space  Administration  TN  D-2181, 
Analytical  Study  of  Aircraft-Developed  Spins  and  Deter- 
mination of  Moments  Reguired  for  Satisfactory  Spin 
Recovery,  by  E.  L.  Anglin,  and  S.  H.  Scher,  1964. 

9.  Champoux ,  R.  L.   The  Numerical  Solution  and  Analysis  of 
Airplane  Spin  Eguations  Modeled  in  an  Inertial  Coordinate 
System,  Masters  Thesis,  United  States  Naval  Postgraduate 
School,  Monterey,  California,  March  1972. 


64 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

2.  U.  S.  Naval  Test  Pilot  School 
Naval  Air  Test  Center 
Patuxent  River,  Maryland  20670 

3.  Flight  Test  Division 
Naval  Air  Test  Center 
Patuxent  River,  Maryland  20670 

4.  National  Aeronautics  and  Space  Administration 
Langley  Research  Center 

Hampton,  Virginia  23365 

5.  Library,  Code  0212 
Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Chairman,  Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

7.  Professor  L.  V.  Schmidt,  Code  57Sx 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93940 

8.  Professor  T.  H.  Gawain,  Code  57 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

9.  LCDR  R.  R.  Buehler,  USN 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

10.   Professor  M.  H.  Redlin,  LT JG ,  USNR 

Department  of  Aeronautics  (Code  57Rp) 
Naval  Postgraduate  School 
Monterey,  California  93940 


65 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 


(Security  classilication  of  title,    body  of  abstract  and  indexing  annotation  must  be   entered  when   the  overall  report  Is   classified) 


1     ORIGINATING    ACTIVITY   (Corporate  author) 

Naval  Postgraduate  School 
Monterey,  California  93940 


2J.   REPORT    SECURITY     CLASSIFICATION 

Unclassified 


2b.    GROUP 


3      REPORT     TITLE 


The  Derivation,  Solution,  and  Analysis  of  Airplane  Spin 
Equations  Modeled  in  an  Inertial  Coordinate  System 


4     DESCRIPTIVE   NOTES  (Type  ol  report  and,inclusive  dates) 

Master's  Thesis;  March  1972 


5     au  T  ho  R  is  i  (Firs  t  name,  middle  initial,  last  name) 


Roy  Robert  Buehler 


6      REPOR  T    D A  TE 


March    1972 


7a.     TOTAL    NO.    OF    PAGES 

67 


76.    NO.    OF    REFS 

9 


Sa.    CONTRACT    OR    GRANT    NO. 


6.    PROJEC  T   NO. 


9a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.   OTHER   REPORT  NOISI  (Any  other  numbers  that  may  be  assigned 
this  report) 


10.    DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited 


II.    SUPPLEMENTARY    NOTES 


12.    SPONSORING    MILITARY    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California  93940 


13.  ABSTR  AC  T 


The  general  equations  of  motion  for  a  rigid  body  are  derived  in 
cylindrical  coordinates  by  Lagrangian  dynamics  and  used  to  model  the 
motion  of  an  airplane  in  a  steady  spin.   After  simplification,  the 
equations  are  cast  into  a  form  utilizing  conventional  aerodynamic 
derivatives  along  with  other  derivatives  which  may  be  significant 
in  spins.   An  iterative  numerical  solution  procedure  is  outlined 
which  should  simplify  the  problem  of  solving  the  nonlinear  differ- 
ential equations,  and  relationships  between  the  Euler  Angles  used  in 
the  equations  and  the  more  familiar  ordered  set  of  pitch,  roll,  and 
yaw  are  derived  to  permit  computer  input  and  output  of  orientation 
to  be  more  easily  visualized.   The  inverse  problem  is  also  considered, 
and  equations  are  derived  to  translate  spin  test  data  into  the  six 
Lagrange  coordinates  used  so  that  computer  calculations  may  be 
checked  against  test  results  and  actual  aerodynamic  forces  can  be 
compared  with  computed  values.   The  paper  concludes  with  a  dis- 
cussion of  the  implications  of  this  spin  analysis  along  with 
applications  and  extensions. 


DD 


FORM 

1    NO  V    05 

S/N    01 01 -807-681 1 


1473 


(PAGE    1 ) 


66 


Security  Classification 


A-31408 


Security  Classification 


KEY     WO  R  OS 


Airplane 

Spin 

Equations 


D    ,F°»R."„1473     'BACK 


'10t -807-6821 


67 


Security  Classification 


A-  3  t  409 


. 


/ 


Thesis 

B8374 

c.l 


133286 

ler 

The  derivation,  solu- 
tion, and  analysis  of 
airplane  spin  equations 
modeled  in  an  inertial 
coordinate  system. 


133286 


Buehler 

The  derivation,  solu- 
tion, and  analysis  of 
airplane  spin  equations 
coordinate  system. 


thesB8374 

T™.,der'Vatl0n'  S0lutl0n-  and  analysis  o 


3  2768  002  07901  4 

DUDLEY  KNOX  LIBRARY 


